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Abstract 

By using properties of the Young orthogonal representation, this paper de- 
rives a simple form for the Fourier transform of permutations acting on the 
homogeneous space of n-dimensional vectors, and shows that the transform 
requires 2n — 2 multiplications and the same number of additions. 
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1. Introduction 

Let Sn denote the symmetric group on n elements, and S*" the subgroup 
fixing the n-th element. This paper derives a simplification for the Fourier 
transform of Sn acting on = {1, 2, . . . , ra}, or equivalently, the coset space 
Sn/ S^. Fourier analysis of permutations on /„ is important for the statistical 
analysis of ranked data pattern matching, and other applications. 

To put the aim of the paper in context, it is useful to consider the ordinary 
Fourier transform. Let Ti be the n x n unitary matrix with entries Ti-k/ = 
(v^)~^e~-'^'^^^^^/". Then X = Tix is the discrete Fourier transform of the 
vector X. If is the translation operator that sends n i— > {n + d) modn, 
and $rf is the phase shift matrix diag[l, e-^^'"''/", . . . , ei27rc((n-i)/nj^ then 



Similarly, the permutation Fourier transform presented below converts per- 
mutations on /„ to group representation "phase" shifts. 
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Fast Fourier transforms on the groups S'„ and their homogeneous spaces 
have been studied previously. In particular, by applying the method of 
Clausen j2[, Maslen and Rockmore 0, Thm 6.5] give an upper bound for 
the number of operations (either multiplications or additions) on Sn/S^ as 
— ri^. Maslen [4|, Thm 3.5] improves the bound on the same space to show 
that, at most, 3n{n — l)/2 operations are necessary. This paper shows that 
2n — 2 operations are sufficient. 

2. Background for this paper 

We use standard results for permutations [sj. An adjacent transposition 
is the permutation = {k,k + 1) that exchanges the k-th and {k + l)-th 
elements but leaves all others unchanged. Every permutation may be written 
as a product of adjacent transpositions. 

The Fourier transform on Sn relies on the group's irreducible unitary 
representations, with "frequencies" given by arithmetic partitions. Let u = 
(ni, . . . , Hq) be a partition of n with > rij+i and rii + . . . + Ug = n; we write 
u \- n. For every u \- n there exists an irreducible representation, denoted 
Dt,. For example, when z/ = (n), we have = 1 for all a G Sn- For 

other z/, we use the Young orthogonal representation (YOR) to construct the 
matrices. The Fourier transform of / : S'„ ^ C is 

F{u)=J2fi^)D.{^), ^^n- (2) 

For each z/, the coefficient F{v) is a rij, x Uy matrix. If /(a) = 
i.e., / and g are left translates of each other, then, in a manner similar 
to ([1]), we obtain that G{u) = D,j{5YF{u). Of particular interest in this 
paper is the "fundamental frequency" of the transform given by the partition 
= (n — 1, 1). The {n — 1)^ entries of are obtained from the YOR as 
described in detail below. 

It suffices to describe on the adjacent transpositions {t/j}, for those 
generate Sn- Let D^{ti) be the (ra— l)-dimensional matrix diag[l, — 1]. 

For any m, let Xm denote the m-dimensional identity matrix, and for k = 
2, ... n — 1, let Rk be the 2x2 symmetric matrix 
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Now, for k = 2, . 
matrix 



. , n — 1, define D^{Tk) to be the symmetric, block- diagonal. 



-^n—k—l 








Rk 

Ik-2 



(4) 



It may be verified that the matrices {D^{Tk)} satisfy the Coexeter relations 
[s], pg 88], and generate the irreducible YOR for partition = (n — 1, 1). 
Furthermore, note that the decomposition of each cr G 5^ into {Tk} excludes 
Tn-i- Therefore, from (jl]), it follows that, with © denoting matrix direct sum 



and 0„^o(a) 



a [n 



2)-dimensional orthogonal matrix. 



,0- 



[a), for a e SI: 
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3. Fourier analysis on the homogeneous space 

Our goal is to simplify ([2]) for functions defined on /„. We may extend 
each / defined on to a corresponding function / on Sn by /(a) = f{a{n)). 
Note that / is constant on left cosets of and, therefore, "band-limited". 

Proposition 3.1. Given any complex-valued function f defined on In, the 
Fourier coefficients {F(z/)} of the function f on Sn defined by f{cr) = f{cx{n)) 
are such that F{v) = unless v = (n) or v = cf) = (n — 1, 1). 

Proof. Since f{a5) = f{a) for 6 G 5^, we have by © that F(z/) = F(z/)D^(5)*. 
By averaging both sides over Sn, we get -F(z^) = -F(z/)Z(z/) where 

zt'') = 7;r^ E DM- (6) 

Now, the Branching Rule [sl, Thm 2.8.3] shows that for z/ = = (n — 1, 1) and 
u = (n), the representation D^, reduces on the subgroup Sn to contain the 
constant representation, and that no other irreducible representation does 
so. By orthogonality, those matrix entries that are not constant on S"" must 
sum to zero over the subgroup. Therefore 2'(z/) = if is not (n) or 0. □ 
If Z{(f))ij is the (i, j)-th element, then from ([5]) we have = 1, and, 

by orthogonality, Z{(j))ij = for all other {i,j). Since F{(j)) = F{(j))Z{(j)) we 
obtain that F{(f)) is zero except possibly in the leftmost column. Hence, the 
Fourier transform ([2]) need only be calculated for the partition (n), and for 
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the n — 1 entries in the left most column of D^. Let JF denote the the linear 
transformation taking any n-dimensional vector x on J„ to its n Fourier 
transform coefficients X{{n)), and the leftmost column entries X{(f))i^i for 
z = 1, 2, . . ., n — 1. We write ^ 

X = J^x (7) 

to express the transform, now viewed a matrix operation. The transform ([7]) 
requires at most ra^ multiplications and n{n — 1) additions. We show below 
that, in fact, 2n — 2 operations of each kind are sufficient. 

Our result relies on the following n x n matrix U, whose shape is similar 
to a "reverse" upper Hessenberg matrix: 



/+1 +1 +1 ■■■ +1 
-1 -1 -1 ■■■ -1 



U 



-1 -1 2 
\-l +1 








+1 \ 

{n-1 



/ 
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Define A = (f/f/*)-^/^, and let 



T = AU. 



(9) 



It is easily seen that T is an orthogonal matrix, and that A is diagonal with 
entries {ak}^=i, with ai = Ai^i = l/y/n, and for > 1, we have 



1 



k,k 



y/{n-k + l){n-k + 2)' 



(10) 



Let X be any complex-valued n x 1 vector, and let X = Tx. To each 
(J E Sni let the n x n matrix P(o") be the permutation matrix obtained from 
the identity X„ with rows permuted by a, i.e., P(cr)jj = Ia{i),j- Note that 
a \—>- P{a) is an antihomomorphism: P{a6) = P{6)P{a). To see that, note 
that for any a we have P{a)ek = ea-i(k) where is [0, . . . , 0, 1, 0, . . . , 0]* 
with 1 in the fc-th position. If x, y are n x 1 vectors, and y = P{a)x, then 
y{i) = x{a{i)) since -P(o") is the permutation operator on column vectors. 
We now establish the following result, comparable to eq. ([1]). 

Theorem 3.2. For every cr G S'„ and all n x 1 vectors x, we have that 



TP{a) x=[l® D^iaY] X 
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Proof. We start by proving for any adjacent transposition that 

TP{n)T' = 1 © D^{Tk) = 1 © D^inY (11) 

Note from IQ, (ITO!) . the m-th row of T for m > 1 sums to zero, with the 
form 

[-am, -Oim, • • • , -ttm, (n - (m - l))am, 0, . . . , 0]. (12) 

The product TPijk) is the same as T but with columns fc, A; + 1 swapped. 
By ( IT2|) . we see that the only rows of TP{Tk) that are affected by the column 
swap are as follows: for = 1, row n is modified; and for A; > 1, rows 
n — {k — 1), n — {k — 2) are modified. Therefore the product TP{Tk)T^ 
is the same as the identity I in all entries with the following exceptions: 
when A; = 1, we have that [TP(ri)T*]^^ = — = —1; and when k > 
1, we have that the 2x2 submatrix, whose upper- left corner indices are 
{n — {k — 1), n — (A; — 1)), has the symmetric form 

-(^ + l)an-(fe-i) (A;2 - l)a„_(fc_i)a„„(fc_2) 1 
_ (A;2 - l)a„_(fc_i)a„_(fc_2) {k - l)al_^^_2) \ 

Subsituting from (|TOl) . we find that the above simplifies to Rk as defined 
earlier in Q, thus verifying ffTTj) for k = 1, 2, . . . , — 1. 

For the general case, note that every cr G may be written as a prod- 
uct of adjacent transpositions cr = r^j ■ ■ -t/j^. Since cr t— > P((t) is an anti- 
homomorphism, we have that 

Pia) = Pin, . . . TkJ = P{nJ ■ ■ ■ P{Tk,). (14) 

Applying a similarity transformation with T yields 

TP{a)r = TP{njr ■ ■ ■ TPinJTK (15) 
On applying (fTTl) we establish the theorem: 

TP(a)r* = [1 © D^iTkJ'] ■ ■ ■ [1 © D^in^Y] = 1 © D^iaY. (16) 

□ 

Note that for the Fourier transform in ([71), we also have 

J^P{a)x = [1 © D^{(r)fX 

from the translation property. Since this is true for all vectors x, we must 
have JF = [AiXi © A2X„_2] T- To see that, note that JF = CT for some matrix 
C, and, by applying the Theorem above, we see that C commutes with all 
matrices 1 © D^; the result now follows from Schur's lemma [5, pg 23]. 
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3.1. Computation of the transform 

The equality T = AU , combined with the matrix structure in ([8]), sim- 
phfies computation. Let ax{n) = x{l), a^in — 1) = x{l) + x{2), . . ., a^il) = 
x{l) + x{2) + - ■ ■ + x{n). Computing all {ax{k)} values requires n — 1 additions 
due to recursion. If X = Ux then X(l) = a2;(l), X(2) = {n — l)*x{n)—axi2), 
. . ., X{n) = x(2) — a^(n). Hence, if has been computed, computing X 
requires {n — 2) multiplies and (^ — 1) additions. Now, since X = T x = AX, 
and A is diagonal, we see that computing X from X requires an additional n 
multiplications. In total, computing the transform X = Tx requires 2n — 2 
multiplications and 2n — 2 additions. Note that computing X = Tx = CAUx 
does not require any extra computation as we may premultiply the diagonal 
matrix C with A. 

4. Conclusions 

This paper describes a simplification of the Fourier transform on Sn/ S^^, 
and shows that the transform requires 2n — 2 multiplications and the same 
number of additions. 
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